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Abstract The crown morphology of Norway spruce ( P/cea ahies Karst) was studied from 617 sample trees. In order to 
model branch and crown descriptors (he linear and nonlinear regression methods were extensively used. Results show that 
the branch length can be (airly well predicted from the distance to the apex of the tree and that the branch spread has a high 
correlation with the branch length and the insertion angle. Models have been set up to predict the crown ratio, height to the 
base of the living crown, height to the first living branch and the height to the first dead branch from the usual whole-tree 
measurements, namely diameter at breast height, total height and total age. 
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Introduction 

Norway spruce (Picea ahies Karst) is one of the most 
important conifer in France. In 1991, the total area of 
Norway spruce was estimated to be 723000 hm 2 . Its 
volume was 127 million m\ The annual volume incre¬ 
ment was 5.74 million nr. It is widely accepted that 
larger initial spacings should be used now. In 1940's, 
the density of the plantation ranged from 5000 to 10000 
stems/hm 2 . Now the initial density varies between 1000 
to 2500 stems /hm 2 . The largest spacing may be 
650-800 stems /Inn 2 . So it has a significance to study 
the branch and crown development of Norway Spruce. 
There are two aims in this study. One is to set up a 
limbsize model that links usual tree measurements to 
the required inputs of a wood quality simulation soft¬ 
ware (SIMQUA: Leban et al, 1990). Another aim is to 
construct models that predict the main characteristics of 
the crown for Norway Spruce. The results of the pres¬ 
ent study will be integrated into a system for predicting 
the quality of the coniferous wood resources from the 
database of the French National Forest Survey (IFN: 
Inventaire Forestier National). 

Materials and Methods 

Data collection 

Three datasets were used. Dataset I contains 141 trees 


of Norway Spruce, four subsamples which have been 
described in Colin and Houllier (op. cit. 1992). Dataset 
2 contains 452 trees sampled and measured by the Insti¬ 
tute for Forestry Development (IDF: Institut pour le 
Developpement Forestier). Dataset 3 contains 24 sam¬ 
ple trees measured in the context of an EEC (European 
Economic Community) contract. Their attributes are 
presented in Table 1. 

Variables 

Three types of variables were used: “branch descrip¬ 
tors", “whole tree descriptors" .and “crown descrip¬ 
tors". The latter two are the usual tree measurements 
and different crown heights and crown ratios. Variables 
are: 

Branch descriptors : 

A=absolute distance from the upper bud scale scars 
of the annual shoot to the top of the stem, cm; 

ANGLE= external insertion angle of the branch with 
the stem, in degree; 

/?/,=branch length, cm; 

/LS'=branch spread, cm. 

Tree descriptors: 

,4=total age of the tree, years; 

D=diameter of the stem at breast height, cm; 

//=total height of the stem, m; 

D/H -ratio between D and //. cin/ciu. 
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Crown descriptors: 

//FC=height to the first contact, m: 

HFLB =height to the first living branch, m; 
HFDB =height to the first dead branch, m; 
//FZ,C=height to the base of the living crown, m; 
CD=niaximum crown diameter, m; 

C/?= 100( H-HFLC )///, %; 

CR 1 = 1 00(H-HFLB)/H, %. 


Table 1. Attributes of 24 sample trees 


age 

D 

11 

CD 

mm 

HFLB 

If FIX' 

HFC 

53 

38.0 

27.25 

4 46 

1 95 

11 95 

1 5.75 

21 75 

53 

38.0 

28 25 

4.56 

2.00 

12.40 

14 10 

19 70 

53 

31.5 

26.67 

3.74 

2 30 

13.85 

16 40 

21.00 

53 

34.0 

27.70 

4.50 

2.15 

15 35 

17.30 

20.55 

53 

26.5 

26.06 

4.32 

1.85 

14 40 

15.5 

23.00 

53 

27.5 

25 48 

3 94 

225 

14 05 

14 90 

21 00 

53 

31.0 

29.13 

4 08 

2.15 

15 30 

19.10 

24.50 

53 

28.5 

28.79 

3.38 

2.05 

19 55 

21 5 

24.55 

53 

30 0 

29 to 

3 94 

2.15 

19 20 

20 15 

22 25 

53 

34.0 

28.56 

3.60 

1.75 

16.50 

19 00 

23 55 

53 

26.5 

28.20 

3.70 

2.30 

15.40 

1935 

25 20 

53 

33.0 

28.84 

5.24 

1.80 

16.20 

I 860 

22 25 

63 

30.0 

25.05 

4,24 

2.40 

10 20 

13 60 

2100 

63 

45.5 

25 24 

4.76 

2.70 

9.05 

II 65 

19 80 

63 

32.5 

25,30 

4 20 

2 15 

9,75 

16,20 

21.75 

63 

35 0 

24.74 

4 64 

2 65 

7.75 

14 35 

18.70 

63 

43,5 

27 80 

4 08 

2.45 

14,80 

16.90 

21 70 

63 

380 

27.30 

4.38 

2 70 

10.20 

16 15 

21 50 

63 

33.0 

23 70 

4.16 

235 

12.45 

14 05 

18 25 

63 

32.0 

26.26 

4.22 

• 2.35 

12 35 

15 55 

20 15 

63 

34 5 

27.07 

3.92 

2.50 

16.05 

17 05 

21.35 

63 

38.0 

26.05 

4.84 

1.90 

13.50 

14.30 

21.05 

63 

31.5 

25.93 

3 70 

2.20 

12.45 

15 30 

21.45 

63 

38.0 

26.07 

4.12 

2.15 

10.85 

15.30 

22.30 


where: C/>=maxiiniim crown diameter, m, 

/7=diameter of the stem at breast height, cm; 

//=total height of the stem. m. 

///■'<-height to the first contact, m; 

///•7,fl=height to the first live branch, m: 

IIF1)B =height to the first dead branch, m; 

///•y.fl=heighl to the base of the live crowti.m 

Models 

In order to select the optimal model, several kinds of 
models were used. Referencing to others studies 
(Mitchell, 1975; Kajihara, 1975, 1976, 1977; Dyeretal r 
1987; Hashimoto, 1990, 1991; Ottorini, 1991; Colin et 
al, 1991, 1992),the following models were used. 


Models for branch length: 

BL-a( 1 -exp( -b/aX)) L (I) 

BL=aX i)J5 -bXexp(-cX) (2) 

BL=a\n(-b/aX+ 1) (3) 

BL=a\ h exp(-cX) (4) 

BL=aexp(-bX) (5) 

Models for branch spread: 

BS=a+bBLsm{AN(]LE) (6) 

BS=(a+bexp(-cBL))BL (7) 

BS=a+b BL (8) 

Models for the height to the base of the living crown 
and to the first living branch: 

HBLC or HFLB=Hexp(-(aD/H+bAH 2 +cH)) (9) 

HBLC or HFLB=Hexp(-(aD/H+bAH 2 +cH+dA )) (10) 

HBLC or HFLB=Hexp(-(aD/H+bAH 2 +cA )) (11) 

HBLC or HFLB=Hexp(-(aD/H+bA ' 5 )H 2 +cA 15 )) (12) 

HBLC or HFLB=H( 1 -7/-exp(-M 15 ) H 2 +cIi/D+dH 2 ) (13) 

Models for the height to the first dead branch : 
HFDB=a+b(H -1 30)/D+cA 1 *H+dA ' 5 Z> (14) 

HFDB=a+bH/D+cA H+dA D (15) 

Models for the crown ratio : 

CR= 100( I -exp (-(aD/H+ bA if+cll))) (16) 

C 'R or CR I = 100(a+exp( -bA 1 -)+cH/D) (17) 

CR or C7? I=I00(a+exp( -bA 1 *)+cII/D+dA 15 D) (18) 
CR\ = 10()( 1 -exp( -(uD/H+hA H 2 +cH+dA))) (19) 


where: a, b, c and d are the parameters to be estimated. 

Statistical analysis 

With SAS (Statistical Analysis System, versions 6.0) 
(SAS Institute, 1988) package, linear and nonlinear 
regression methods have been extensively used. The 
proposed equations were chosen as compromises be¬ 
tween (i) the search of a good fit as measured by ad¬ 
justment statistics and by a visual analysis of residuals 
and (ii) the parsimony and the robustness of the model 
(i.e. we tried to avoid a too great number of parameters). 

The following results include parameter estimates, 
their standard error, and their 95%-confidence interval, 
root mean squared error (RMSE) or weighted mean 
squared error (WMSE). adjusted Rsquare (R^- 1 -((»- 
1 )/(n-p))( 1 -R 2 )) and weighting expressions. 

Results 

Crown shape and size relations 

Crown shape and size result from the relationship be¬ 
tween crown extension (i.e. branch growth) and height 
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growth. The branch growth involves in branch length 
(/>’/.).branch spread) BS) and branch diameter. Branch 
diameter is not discussed here. 

Branch length ( BL ) 

By titling equations (1) to (5), we have found that equa¬ 
tion (I) is the best n tree top to the first contact, to the 
base of the living crown and to the first presented in 


Table 2. 

Branch spread (BS) 

Branch spread has a linear relationship with branch 
length and the sine of the insertion angle. By computing, 
equation (6) is the best. Parameter estimates for the 
equation (6) are given in fable 3. 


Table 2. Parameter estimates for models of branch length 


Position 

Dataset 

n 

Weight 

Parameter 

Asymptotic 

Asymptotic confidence interval 95% 

WMSF 

i, 2 





estimates 

standard ciror 

lower 

upper 







<>-=2X7 4455 

6X154 

274 03X3 

300 8508 




dataset 3 

340 

\l.x" 

btl 4403 

0 WO 

0 1554 

1 7251 

1 252 

0 93 

Inst cimtuct 




c= 1 4500 

0 6224 

0.2257 

2 6743 



eq ( 1 ) 




</-296 X7I2 

17 455.3 

262 6250 

.3.31 1175 




dataset 3 + dataset 1 

1220 

I/.V 1 ' 

b-U 7334 

0.0X55 

0 5656 

0 9012 

1.144 

0 83 





c~~ 1 1W5 

0 ! 159 

0 962! 

1 4170 







(!--• 312 6960 

8 424! 

296 1497 

329 2423 




dataset 3 

574 

Ur 1 * 

h ~0 9027 

D0350 

0 8339 

0 9716 

1 203 

0 95 

bnse of crow n 




c~ 1 3998 

0 0440 

1 3135 

1 4862 



cq(2> 




o'271 I2bb 

5 5640 

260 2138 

282 0.394 




dataset 3 1 dalasel 1 

2406 

l/x 

b- l) 8104 

0 0322 

0 7472 

(! 8736 

0 905 

0 84 





rl 2739 

0 04.3.5 

1 18XX 

1 3592 







f r.303 (,293 

5 55111 

292 7320 

314 5266 




dataset 3 

w 

l/.f' 

/>“ () Q* | (a 

1)0.317 

0 8713 

0 9918 

1 033 

0.95 

fust living branch 




c ~I 4339 

(70400 

1 3555 

1 5124 



cq 1 2 ) 




0-272 0651 

3 7259 

264 7594 

279.3708 




dataset 3 + dataset 1 

290b 

1 lx 

/>—(I.RI06 

0 0261 

0.7594 

0.8619 

1.089 

0.85 





i-l 2744 

0 0369 

1.2024 

1 3592 




from a static point of view, equations (I) and (6) arc growth , these equations reflect the process of expan- 

an expression of crown shape and size. As for a give A sion of the parts of a crown, 

varies with tree height in association with height 


Table 3. Parameter estimates for the prediction of branch spread (equation (6)) 


Dataset 

Weight 

Parameter 

Asymptotic standard error 

WMSF. 

^adj 



estimates 




dalasel 3 

l/(/;/.sin(.rV(,7/D) 

<i~0 81528 

000393 

6 63(1 

0.98 

dataset 1 

l/<«/.sin(.1,Y<;/./•.)) 

«-() 75614 

0 00465 

17 376 

0.99 

dataset 3 t dataset 1 

l/t«/.sin(,l,V<;/.A)) 

<r7 8821 

0 00324 

4.719 

0.98 


Models of main crown characteristics 

Global description of the crown The depend¬ 
ent variables were height to the lirst dead branch 
(IlFDB ).height to the first living branch (IIFLB). height 
to base of the living crown (HBLC) and crown ratio 

(c/e). 

The independent variables were total height (//).total 


age (/I (.diameter at breast height (D) and various com¬ 
binations of these variables such as: 1///, IF. H/D , D/H. 
A II. A I). etc.... 

Height to the base of living crown (HBLC) The 
best individual predictors are A, D/H and //. By fitting 
equations (9) to (13), equation (9) is the best model for 
HBLC. In order to take into account the fact that the 
data set includes both data for isolated trees and data 
for trees belonging to die same stand, the weight of 


to 
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each tree was inversely proportional to the number of 
trees belonging to the same stand. This weighting pro¬ 
cedure led to a good fit especially for the data collected 
on old isolated trees.Parameter estimates are presented 
in Table 4. 

Height to the first living branch(//FLB) We 


used the same method (equation and weighting expres¬ 
sion) as for HBLC. We have found equation (10) is the 
best model for predicting HFLB. As previously, the 
weight expression took into account the number oi 
sampled trees in each stand. Parameter estimates are 
showed in Table 5. 


Table 4. Parameter estimates for the prediction of height to the base of the living crown (equation (9)) 


Dataset Parameter estimates Asymptotic standard error _ Asymptotic confidence interval ^5% _ 

__lower_upper_ VVMS1- _ ^adj 


((-68 90176 

5. #4732 

57 33747 

80 46606 



dataset 1 6=-5 6* 10"’ 

6 2* 10 " 

-6 8*! O'" 

-4 4* 10" 

0 (MM) 

0 75 

C--O.OOI47 

0.00397 

-0 00932 

0 00638 



o-35 61325 

4 10335 

27 54785 

4367865 



dataset 2 />=-6 0*It) ' 

9 0*! ()" 

-7 7* nr' 

-4.2*10'* 

0,001 1 

0 61 

c=0.04114 

0 00677 

0 02862 

005366 



0=62,20247 

2 4451 | 

57 40020 

67 00474 



all* 6=-5 8* 10" 

3.5*1 O' 7 

-6 4*1(4" 

-5 | * I(> '• 

0 0030 

0 79 

(=0.00209 

0 00180 

-0.00145 

0 00563 



* all-dataset 1 <■ dataset 2 *■ dataset 3 


Table 5. Parameter estimates for the prediction of height to the first living branch (equation (10)) 


Datase Parameter estimates Asymptotic standard error 

Asymptotic confidence interval 95% 





lower 

upper 

wmsp: r«i, ! 


o=IOI 63831 

II 12893 

79 62850 

123 6481 



dataset 1 />=-8 s*|()= 

2 7* Iti 

-1 3*10"' 

-3 2* 1 O ' - 



(=0 00356 

0 00672 

-0 01684 

000973 

0 0106 0 58 


(7=0 1)0296 

0.00227 

-OOOI54 

0 00745 



(/=|04 38967 

10 20408 

84 23645 

24 5429(1 



6=-4 l*IOY, 

2 4*1 O' 7 

-1 4*10" 

-4 3*10" 



dataset 1 + dataset 3 ( =000683 

0 00594 

-0.01856 

0 00490 

00094 0 61 


rf-0 00371 

0.00210 

000044 

000787 



Height to the first dead hranch(/yFWB) The 

statistical analysis was carried out on 80 trees in data 1 
(primed trees were removed). Models for HBLC and 
HFLB were first tested but the best results were ob¬ 
tained with a linear model (equation 14) including H A, 
H/D and D A. Parameter estimates for equation (14) are 
given in Table 6. 

Crown ratio(C7?) Since CR= 1 Q()(H-HBI.(')/!/ 

and CR I = 1 ()0 ( H-HFLB)/H and best models for HBLC 

and HFLB are equation (9) and equation (10) respec¬ 
tively. we finally got equation (16) and equation (19). 
By fitting equations (16 to (18), equation (16) is the 
best model for CR. By fitting equations (17) to (19), 
equation (19) is the best model for CR\. Parameter 
estimates for equation(l6) and equation (19) are the 
same as in table 4 and in Table 5 respectively. 

Table 6. Parameter estimates for the prediction of height to the first dead branch (equation (14)) 



Parameter estimates 

Asymptotic standard error 

Asymptotic confidence interval 95% 



Dataset 


lower 

upper 

wmsp: K ad| 2 


a--11 58979 

36 65005 

-150,58491 

-4 59466 



b 1 64623 

dataset 1 (--5 7*1 o ' 7 

0 47344 

3 2*10 ^ 

o 70.328 

-4 

-1 2*10 

2 58918 

-6 

7.3*10 

0.00081 0 50 


(/= 0.00023 

0 00006 

-0 00011 

0.00036 
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Conclusion 

These models provide a comprehensive and consistent 
framework for predicting branchiness and crown char¬ 
acteristics from usual tree measurement as those avail¬ 
able from local management or the database of the 
French National Forest Survey. 

Above results show that: (1) The branch length has a 
high correlation with the distance from the tree top. The 
size of whorl branch may be predicted at any point 
along the stem; (2) The main characteristics of the 
crown may be estimated from usual tree measurements, 
such as total age. total height and diameter at breast 
height, etc.. (3) The weightiness is available on model¬ 
ing the shape of the crown. 

References 

Baxter. S.M. & Canned. M.G.R. 1979. Branch development on 
leaders of I’icea sitchcnsis. Can. J. For. Res.. 8:121-128 
Beadle. CM,.. Tablot. II. & Jarvis. P.G. 1982. Canopy structure 
and leaf area index in a mature Scots pine forest. Forestry. 
55:105-123 

Bouchon, J. 1984. Importance des plantations de Douglas et 
Epic a cn France. Revue Forcsticre Fran aisc. 26. 254-258 
Colin. F. & I loullicr. F. 1991. Branchiness of Norway spruce in 
north-eastern France: modeling vertical trends in maximum 
nodal branch size. Ann. Sci. For.. 6:679-693 
Colin. F. & I loullicr. F. 1992, Branchiness of Norway spruce in 
north-eastern France: predicting the main crown characteris¬ 
tics from usual tree measurement. Soumis mix Ann. Sci. For. 
Curtis. R.O. & Reukema l).l„ 1970. Crown development and site 
estimates in a Douglas-IJr plantation spacing test. Forest Sci¬ 
ence. 16(3 ):287-30l 

Dyer. M l .. & Burkhardt III. 1987. Compatible ratio and crown 
height models. Can. J. For. Res.: 17:572-574 
Ilashimoto. R. 1990. Analysis of the morphology and structure 
of crowns in a young sttgi ( Cryptomeria japonica ) stand 
tree physiology. 6:119-134 

Ilashimoto. R. 1991. Canopy development in young sugi 
( Cryptomeria japonica ) stands in relation to changes with 


age in crown morphology and structure. Tree physiology. 
8:129-143 

I loullicr F„ Bouchon. .1. et Birot. Y. 1991. Mod lisation de la 
dynamique dcs peuplemcnts forestiers: 1 tat et perspectives. 
Revue Forestiere Fran aise. XLIII(2):87-108 
Kajihara. M. 1975. Studies on the morphology and dimensions 
of tree crowns in even-aged stand of sugi, J. Crown form ( in 
Japanese with English summary ). .1. Jpn. For. Soc.. 57:425- 
431 

Kajihara. M. 1976. Studies on the morphology and dimensions 
of tree crowns in even-aged stand of sugi. 111. Development of 
crown morphology with growing stage ( in Japanese with 
English summary ). .1. Jpn. For. Soc.. 58:313-320 
Kajihara. M. 1977. Studies on the morphology and dimensions 
of tree crowns in even-aged stand of sugi. V. Crown surface 
area and crown volume ( in Japanese with English summary ), 
.1. Jpn. For. Soc. 59:233-240 

Ecban. J.M. el Duchanois G. 1990. SIMQUA: Un logiciel de 
simulation de la qualit des bois, Ann. Sci. For. 47(5). 483- 
493 

Madgwick llai, Tann Co et Fu Mao-Yi, 1986. Crown develop¬ 
ment in young I’icca abies stand. J. For. Res.. 1:195-204 
MAGUIRE. D.A. & I1ANN. D.W. 1987 A stem dissection tech¬ 
nique for dating branch mortality and reconstructing past 
crown recession. Forest Science. 33(4):858-87l 
Mitchell. K..I. 1975. Dynamics and simulated yield of Douglas 
Hr. Forest Science. Monograph 17. 39p 
Ottorini. J.M. 1991. Growth and development of individual 
Douglas-fir in stands for applications to simulation in silvicul¬ 
ture. Ann. Sci. For.. 48(6):651-666 
SAS Institute. 1988. SAS/STAT. Guide for personal computers. 

1028p. 

Satoo. T. and Imolo II. 1979. Modelling crown canopy of an 
even-aged stand of Cryptomeria japonica from measurement 
of leaf mass: A new approach to the ntophology of forest 
crown. .1. Jpn. For. Soc. 61:127-134 
Tome. M. & Burkhart. II. IE 1989. Distance-dependent competi¬ 
tion measures for predicting growth of individual tree. For. 
Sci. .35: 816-831 

(Responsible Editor: Chai Rnihai) 


12 



